Surface entanglement in quantum spin networks 
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We study the ground-state entanglement in systems of spins forming the boundary of a quantum spin network 
in arbitrary geometries and dimensionality. We show that as long as they are weakly coupled to the bulk of 
the network, the surface spins are strongly entangled, even when distant and non directly interacting, thereby 
generalizing the phenomenon of long-distance entanglement occurring in quantum spin chains. Depending on 
the structure of the couplings between surface and bulk spins, we discuss in detail how the patterns of surface 
entanglement can range from multi-pair bipartite to fully multipartite. In the context of quantum information 
and communication, these results find immediate application to the implementation of quantum routers, that is 
devices able to distribute quantum correlations on demand among multiple network nodes. 

PACS numbers: 03.67.Bg,75.10.Jm,03.67.Hk 



I. INTRODUCTION 

Quantum spin networks are being actively investigated for 
quantum technological applications ranging from quantum 
teleportation to the quantum internet. Their dynamical prop- 
erties have been extensively studied in order to assess their 
functionality for quantum state engineering and quantum state 
transfer IUl4l5ll. while their static properties have been ana- 
lyzed in order to identify ground or equilibrium thermal states 
useful as entangled quantum resources lfl6l - l25ll . In principle, 
they can be simulated with a va riety of atom-optical, molec- 
ular, and solid-state platforms ||26l - l32ll . which makes them 
extremely flexible models for the implementation of quan- 
tum technologies. On the other hand, a very important task 
in many fundamental protocols of quantum information and 
communication is the harnessing of entanglement between re- 
mote quantum objects. In particular, the distribution of en- 
tanglement between distant, not directly interacting systems, 
which may be part of a network, have bee n propo sed and ex- 
plored in many different physical settings Il2ll425l 133144611 . In 
this context, a relevant concept is that of long-distance entan- 
glement (LDE) I12ll - l25n . namely the emergence of a sizeable 
end-to-end entanglement in spin chains whenever the cou- 
pling of the end spins to the remainder of the system is suffi- 
ciently weak. The occurrence of this intriguing phenomenon 
in one-dimensional systems naturally suggests to investigate 
the existence of analogous or more general effects in higher 
dimensions and nontrivial geometries. 

In the present work we extend the notion of LDE 1I21I - I25I1 to 
higher dimensional lattices and generic geometries. We intro- 
duce the concept of surface entanglement, which refers to the 
emergence of entanglement in the reduced ground-state den- 
sity matrix of the boundary spins of the network, that is the 
surface spins. In this respect, the LDE can be seen as a par- 
ticular case of this more general phenomenon, corresponding 
to the case in which the surface coincides with two end spins 
of a liner spin chain. As we shall see, the main feature of 
surface entanglement, at odds with LDE, is that it can involve 
a large number of distant, non-interacting spins so that many 
different structures of the shared correlations can be identi- 
fied between them, ranging from large collections of bipartite- 
entangled spin pairs to complex patterns of strong multipartite 



entanglement. On the other hand, similarly to LDE, one of 
the conditions for the onset of surface entanglement is that the 
couplings between the spins of the surface and the ones in the 
bulk be much smaller than the coupling strengths among the 
spins in the bulk. This feature implies that surface entangle- 
ment can be controlled by adjusting locally the interactions 
between the surface spins and the bulk, while the actual inter- 
actions within the rest of the network are up to a certain de- 
gree (as suitably specified below) irrelevant. This makes the 
present approach appealing for the realistic implementation of 
quantum routers llM l47l - l5lll . namely devices able to distribute 
quantum correlations between distant nodes of a network of 
quantum systems. 

The paper is organized as follows. In Sec. [II] we introduce 
the class of quantum spin models to investigate, and we derive 
an effective Hamiltonian for the surface spin-dynamics in the 
limit of small coupling between surface and bulk spins, show- 
ing that in this limit the effective Hamiltonians share always 
the same symmetries of the original models. In Sec. IIII Al we 
address the most elementary case in which the surface is made 
of only two spins, an immediate generalization of the LDE in 
linear chains to bulks of arbitrary dimensionality and geome- 
try. In Sec. IIII B I we consider the general situations in which 
the surface is made of an arbitrary number of spins, and we 
discuss how the tailoring of the interactions between surface 
and bulk spins can yield either large arrangements of bipar- 
tite entangled spin pairs on the surface or complex patterns 
of multipartite entanglement involving many or even all of the 
surface spins. Finally, in Sec.|IV]we draw our conclusions and 
discuss possible outlooks. 



II. SURFACE SPINS: TOTAL AND EFFECTIVE 
HAMILTONIANS, AND SYMMETRIES 

In this paper we analyze the properties of ground-state en- 
tanglement in various classes of spin- 1/2 models for lattices 
of generic geometries and dimensions. We consider structures 
that can be characterized by two subsets of spins, respectively 
B and S , where B indicates the set of bulk spins and S is the 
set of surface spins. These sets are not arbitrary; the set B 
of bulk spins is assumed to be described by an Hamiltonian 
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whose ground state is non-degenerate, thus implying that the 
total number of spins in the bulk is even; the set S of surface 
spins is assumed to be such that the spins in S do not interact 
among themselves, and each spin in S interacts with a single 
spin in B with a strength that is significantly smaller than the 
typical one between the spins in B (see Figs.Q]and[3]for exam- 
ples). In order to be as general as possible we assume that the 
spins are coupled by anisotropic Heisenberg XYZ exchange 
interactions. This ample class of models includes, as particu- 
lar cases, the Ising, XY, XX, XXZ, and isotropic Heisenberg 
XXX models. Therefore, all the systems which we will dis- 
cuss in the following are described by a total Hamiltonian of 
the following form 



Ht — Hr + Hs 



(1) 



where Hb accounts for the interactions between the spins in 
the bulk, and H$b describes the interactions of the surface 
spins with the spins in the bulk. They are defined as 



^ = ZZ 4 <A 

j.keB ae{x,y,z\ 
jeS ae{x,y,z] 



(2) 
(3) 



where <x£ (cr^ ) with a e {x,y, z} are the standard spin opera- 
tors for the surface (bulk) spins. All the coupling strengths J" k 
and KJ are of the same order of magnitude. On the other hand, 
the bulk-surface interactions are weighted by the coefficients 
Aj that can take any positive value in the range [0, 1]. 

An important preliminary result is that in the limit of Aj <sc 
1, the dynamics of the surface spins can be approximated by 
an effective interaction Hamiltonian between them that shares 
the same symmetries (XY, XYZ, XXZ, or XXX) of the original 
model Eq. (HJ. In order to prove this property, let us consider a 
model for which all the weighting coefficients coincide: Aj = 
A Vy e S , and let us define the set {\tpe)} of eigenstates of Hb 
and the corresponding eigenvalues E{ (which satisfy Hb\4>c) = 
Ee\(pe)). Under the assumptions that the ground state |0o> of 
Hb is non-degenerate and that the gap between the ground 
and first excited state is much larger than AKJ one can resort 
to perturbation theory in the parameter A in order to study the 
low-energy eigenstates of the total Hamiltonian Ht = Hb + 
Hsb- At zeroth order, the ground space of Ht is degenerate 
with degeneracy equal to the dimension of the Hilbert space 
for the surface spins, and it is spanned by the states of the form 
\<Po) ® Wp), where the states \i/r p ) form a basis in the Hilbert 
space of the surface spins. Thereby, the effective interaction 
Hamiltonian at second order in A for the surface spins reads 
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Taking into account that the anisotropic Heisenberg Hamil- 
tonian Hb in the bulk preserves the parity with respect to 
the three fundamental directions, i.e. [H B ,P a ] = with 
P a - WjO^ an d a e y> z}, we can rewrite H e g as 



ae{x,y,z] j,keS 



(5) 



where the effective coupling strength h" k is equal to 

A" = -2A 2 K"K"Y ^ . (6) 

We see that this effective Hamiltonian describes, in general, a 
fully connected XYZ model (or XY model if K\ = V£) for 
the surface spins, that is each surface spin interacts with all 
other surface spins. 

Moreover, if the bulk Hamiltonian Hb commutes also with 
the total magnetization along the z-axes, S z = Yjj^B' ^ or 
instance if it is a Hamiltonian of the XXZ form, then the coef- 
ficients A v , and A\, can be rewritten as 

],k ],k 



._ x ^o|o-+ \(f> e ){(f>e\o- „ |0o> 
A* = -2A 2 K]Kl V ^ , 



1*0 
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A ^ = - 2A /^Z ^ — — • (7) 



E{ - Eq 



where the <x+ and cr Bt are the raising and lowering spin op- 
erators, and A > , = t^A*,. Thus, if also Hsb is of the XXZ 

J,K ^j^k J> 

type, with Kl = K?Vk, then also the effective Hamiltonian 
H e ff is of the XXZ type. Finally, in the case that the original 
total Hamiltonian Ht belongs to the class of isotropic Heisen- 
berg XXX models, which commute with the total magnetiza- 
tion with respect to any axes, applying the same reasoning as 
in the previous cases, it is fairly straightforward to show that 
also in this instance the effective Hamiltonian is of the XXX 
type. 

In conclusion, the effective interaction Hamiltonian H e g for 
the surface spins shares the same symmetries of the original 
total Hamiltonian Ht, and it is in general fully connected, re- 
gardless of the range of the interactions in the original model. 
As a consequence, on the one hand, the ground-state prop- 
erties of H e g will be dictated by the symmetries of Ht and, 
on the other hand, given the fully connected nature of H e g, 
one can anticipate that in general its ground state will exhibit 
complex patterns of bipartite and multipartite entanglement. 
These qualitative picture is confirmed quantitatively by the re- 
sults of an extended numerical analysis performed directly on 
the original model Ht and reported in Section Hill 



HI. NUMERICAL RESULTS: EXACT DIAGONALIZATION 
OF THE TOTAL HAMILTONIAN 

In this section we analyze the entanglement properties of 
the reduced ground-state density matrix of the surface spins, 
obtained by tracing out the degrees of freedom of the bulk 
spins after exact numerical diagonalization of the total Hamil- 
tonian Ht for different dimensions, geometries, and sym- 
metries. We will also analyze in detail the energy gap be- 
tween ground and first excited state, because this quantity is 
of central importance in assessing the thermal stability of the 
ground-state entanglement, in particular of that of the surface 
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FIG. 1: Top: two possible simple schemes of spin networks with 
a surface consisting of only two spins. Left: a two-dimensional 
square configuration. Right: a three-dimensional cubic configura- 
tion. Spins are denoted by circles, while the lines connecting the 
circles indicate the spins that interact directly. The bulk spins in- 
teract with a rescaled dimensionless coupling normalized to unity. 
The couplings A identify the two surface spins (5 i,52), i.e. the 
(not directly interacting) spins at the boundary of the network that 
interact weakly with the bulk. Bottom: Concurrence between the 
surface-spins 5 1 and 52 as a function of A for the square (black lines) 
and the cubic (gray lines) configuration. These results correspond 
to XX interactions (solid lines) with J x k = J y k = K x ( = K y ( = J, 
Vy, k 6 B and Vf 6 5, and XXZ interactions (dashed lines), with 
2J\ = 2K) = P = f jk = K} = K] = J, Vj,k e B and W 6 5. 
In both cases the interactions are antiferromagnetic (J > 0). Inset: 
energy gap as a function of the two-spin surface concurrence. 

spins. Furthermore, in the perspective of an experimental re- 
alization of surface entanglement based on adiabatic ground- 
state preparation protocols [52], the size of the gap sets a limit 
to the maximum allowed velocity for adiabatic manipulation. 



A. Two-spin surface 

We consider first the simplest situation in which the surface 
is made of only two spins. This situation can be seen as a 
generalization of the LDE in quantum spin chains 1I2II - I24I1 to 
networks whose bulk has an arbitrary geometry but a minimal 
number (j = 1 , 2) of surface spins, defined as the - non directly 
interacting - boundary spins that interact weakly with the bulk 
via couplings Aj. In the limit of Aj «: 1, and in agreement 
with the findings of the previous section, the two-spin ground- 
state reduced density matrix approaches a maximally entan- 
gled Bell state, with the exception of the cases in which the 
ground state of the two-spin effective Hamiltonian is degener- 
ate, as in the Ising and the ferromagnetic isotropic Heisenberg 
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FIG. 2: Concurrence in the reduced ground state of the surface spins 
5 1 and 52 for square (black line) and cubic (gray line) network con- 
figurations of Fig. [TJ with XXZ interactions (Jj k = J y k = K x ( = K y = 
J > 0, V j, k e B and W 6 5 ), as a function of K z = K' ( ,\/{€S. Both 
curves are plotted for A = 0.1. 



models. Examples of results for the two-spin surface entan- 
glement are reported, for two- and three-dimensional systems, 
and for XX and XXZ interactions, in Fig.Q] We observe that in 
all cases, when Ai = A% = A is sufficiently small, the concur- 
rence of the reduced density matrix of the two surface spins 
approaches unity, meaning that in this limit the reduced state 
correctly approaches a Bell state as was to be expected from 
the discussion in the previous section. Moreover, since the 
effective coupling strengths A" 2 are proportional to A 2 , the 
energy gap tends to vanish as A is reduced, thus making the 
surface entanglement extremely unstable against thermal fluc- 
tuations in the limit of vanishing A. On the other hand, as the 
concurrence decreases for larger A, the energy gap increases 
correspondingly. 

One notices from Fig. Q]that, depending on the geometry, 
an anisotropy in the spin-spin interaction along the transverse 
z direction may be instrumental to the surface entanglement. 
In particular, the cubic network configuration with XXZ in- 
teractions possesses larger entanglement with respect to the 
same geometry with XX interactions. In general, the optimal 
value of K z which maximize the entanglement is geometry- 
dependent. For example, as shown in Fig. [2] at fixed A and 
in the case of the cubic configuration of Fig. [1] the concur- 
rence reaches its maximum, at K z ^ J, i.e. for XXX isotropic 
Heisenberg interactions, whereas for the square configuration 
it is maximal at K- ~ 0.065 /. 

For assigned geometry, symmetries of the interactions, and 
surface-to-bulk coupling A, the concurrence slowly decays 
with the dimension of the bulk (i.e. the number of spins in 
the lattice, excluding the two surface spins). This is shown 
in Fig. [3] where the surface concurrence is reported for di- 
verse models whose bulks have a ring-like shape with increas- 
ing number of spins. Maximum entanglement is obtained 
for the smallest (in terms of number of spins) configurations. 
Viceversa, the gap mildly increases with the number of spins. 
Therefore, this geometry results promising for an actual im- 
plementation of a rooter based on surface entanglement: de- 
pending on the desired working point, it entails a constrained 
optimization of the two contrasting requests of significantly 
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FIG. 3: Top: ring-shaped quantum spin network of increasing size. 
Bottom: surface concurrence (upper panel) and energy gap (lower 
panel) as functions of the total number of spins for this type of net- 
work. The parameter A denotes the weak coupling between the two 
surface spins and the bulk of the network. The interactions are of XX 
ferromagnetic type with J*, k = J y . k = = tC ( = J,Vj,ke B and 
It < S . 

sizeable entanglement and energy gap. 

It should be noticed that the models studied in Figs. Q] 
and[3]satisfy the generalized Toulouse criteria for frustration- 
free systems, as introduced in ll53ll . and therefore they are 
not geometrically frustrated. In general, a model that does 
not satisfy the generalized Toulouse criteria, and that there- 
fore is geometrically frustrated, has a surface-spin concur- 
rence which decays much faster with increasing surface-to- 
bulk coupling A, as compared to the corresponding non frus- 
trated model. This behavior is illustrated in Fig.|4]where mod- 
els that do not satisfy the generalized Toulouse criteria are 
realized with anti-ferromagnetic nearest-neighbor interactions 
on square and pentagonal bulk geometries with an odd num- 
ber (five) of interaction bonds. These are compared with the 
corresponding ferromagnetic models, which satisfy the gen- 
eralized Toulouse criteria and are thus geometrically unfrus- 
trated [53]. Therefore, quantum spin networks with geometri- 
cally non frustrated bulk configurations should always be pre- 
ferred for an optimization of the corresponding surface entan- 
glement. 

Finally, before considering the general case of networks 
with an arbitrary number of boundary spins, we turn to the 
concept of modular entanglement recently introduced for spin 
chains [25], and we look for generalizations to higher dimen- 
sions and generic geometries. It has been demonstrated in 
Ref. B25I1 that, in the case of linear spin chains, the properties 
of the end-to-end entanglement, including its stability against 
thermal fluctuations, can be greatly enhanced when various 
sub-chains are serially coupled by weak interactions to form 
a modular structure. Moving to spin networks with only 
two surface spins, a similar behavior is observed in Fig. 



FIG. 4: Top: examples of two-dimensional networks with bulk con- 
figurations allowing for geometrical frustration. Left: square with 
five interaction bonds. Right: pentagonal with five interaction bonds. 
Bottom: concurrence between the surface spins 5 1 and S 2 as a func- 
tion of the surface-to-bulk coupling A. The interactions are of the XX 
type (/* = J y Jk = K* t = K] = J,Vj,ke B and W e S). Solid lines: 
antiferromagnetic, geometrically frustrated case. Dashed lines: fer- 
romagnetic, geometrically unfrustrated case. Black lines: network 
with square geometry of the bulk. Gray lines: network with pen- 
tagonal geometry of the bulk. Inset: corresponding behavior of the 
energy gap as a function of A. 



where we consider an example in which the bulk of the net- 
work is organized in linearly coupled two-dimensional mod- 
ules (blocks). Indeed, the associated surface entanglement is 
enhanced as the number of modules is increased, while the 
energy gap correspondingly decreases. The mechanism of 
modular entanglement can thus be effectively generalized to 
higher-dimensional geometries. Comparing Fig. |5]and Fig. [3] 
one sees that ring-like and modular networks exhibit oppo- 
site behaviors of surface entanglement and energy gap as a 
function of the number of bulk spins (bulk modules). This 
observation suggests the alternative use or the combination of 
the two configurations depending on the assigned physical or 
communication task. 



B. Many-spin surfaces 

So far we have considered networks whose surfaces con- 
sist of only two spins. When the boundary contains more than 
two spins, then many different patterns of ground-state surface 
entanglement can be identified. In particular, in the follow- 
ing we will consider two relevant limiting cases. Firstly we 
shall study the situation in which the surface spins dimerize 
in a series of strongly entangled pairs (bipartite entanglement 
replicated in a large number of entangled pairs). Secondly, 
we address the case in which the boundary is characterized 
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FIG. 5: Top: Quantum spin network with a modular structure of 
coupled two-dimensional blocks of spins. Bottom: for this type 
of modular network, ground-state concurrence between the surface 
spins 5 1 and 5 1 (upper panel) and corresponding energy gap (lower 
panel) as functions of the total number of blocks for different values 
of the weak coupling A. The interactions are of the XX type with 
/* = J y jk = K* = K y ( = /, V j, k e B and V£ e S , and the same results 
hold for ferromagnetic and antiferromagnetic interactions. 

by strong multipartite entanglement involving essentially all 
of the surface spins. These two limiting cases are important 
because most of the intermediate possibilities, where complex 
patterns of bipartite and multipartite entanglement may coex- 
ist simultaneously, exhibit essentially various aspects of these 
two instances. In order to clarify this complex phenomenol- 
ogy we discuss separately the two limiting cases. 

1. Many-pair bipartite surface entanglement 

Many entangled spin pairs are realized in the reduced 
ground state of the surface spins when their interaction with 
the bulk follows a two by two pattern such that every two 
surface spins share the same coupling A which in turn dif- 
fers from that of all other pairs, as illustrated in Figs. |6]and 
Q This result can be understood as an extension of the situa- 
tion holding for two-spin surfaces - discussed in Sec. MI Al - to 
a series of many nested two-spin surfaces where each added 
term is coupled to the bulk with decreasing coupling strength. 
Let us, for example, consider the case of four surface spins, 
such that two of them are coupled to the bulk with a strength 
A, whereas the coupling strength of the other two is A', and 
they satisfy the relation X «: A <K 1. In complete agree- 
ment with what we have proved in Sec. MI Al the two spins 
with the weakest interactions A' become strongly entangled. 
Indeed, they can be thought as forming the surface of an en- 
larged bulk consisting of the actual bulk and the two remain- 
ing spins in the surface. In turn, also the enlarged bulk can 
be seen as a reduced system with two surface spins that inter- 




FIG. 6: Top: single-square (left) and nested tow-squares (right) two- 
dimensional architectures for quantum spin networks. Bottom: con- 
currence in the reduced ground state of the pairs S1—S2 (solid lines) 
and 53-54 (dashed lines) of surface spins for the single-square 
(black lines) and the nested two-squares (gray lines) configurations. 
The interactions are of the XX type with = J y . k = KJ = K y ( = J, 
~ij,k 6 B and W 6 5 , with the same results holding for ferromag- 
netic and antiferromagnetic couplings. The surface spins 5i and 52 
are coupled to the bulk with strength A, while the surface spins 53 
and 54 with strength A' = A 2 . Inset: behavior of the energy gap as a 
function of A. 



act with the bulk with strength A, and thus also the reduced 
ground state of this pair is strongly entangled. In conclusion, 
the total surface of the network separates (dimerises) in two 
entangled spin pairs that become maximally entangled in the 
limit of vanishing surface-to-bulk couplings. This result ex- 
tends straightforwardly to surfaces with many spins, resulting 
in a structure of many entangled dimers with the inclusion of 
an arbitrary number of spin pairs in the surface of the network. 

We have verified this picture thoroughly by exact numerical 
diagonalization and the results are summarized in Figs.|6]and 
[7] The results for the basic configuration with four surface 
spins are shown in Fig. [6]for interactions of the XX type. The 
concurrence of the surface spin pairs is reported as a function 
of the two weak couplings A for the first pair and A' < A for 
the second pair. Both pairs approach a maximally entangled 
Bell state in the limit of a vanishing coupling strength, thus 
demonstrating the potentiality of the surface entanglement as 
a resource for the parallel distribution of quantum correlations 
among many nodes of spin networks. The multi -pair bipartite 
entanglement occurs also when several pairs are considered, 
as illustrated in Fig. [7] where we consider up to four entan- 
gled pairs in the reduced ground state of the network's sur- 
face, and we report the bipartite entanglement of four differ- 
ent pairs of surface spins with couplings in decreasing order 
A > A' > A" > A'". We see that as the values of the cou- 
plings are progressively reduced in hierarchical order, all the 
two-spin concurrences approach unity and tend to realize a 
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FIG. 7: Top: Ring-shaped network. Bottom (upper panel): for this 
type of network, the concurrences in the reduced ground states of 
the different surface spin pairs: Si — S2 (circles), S3 - S4 (squares), 
55 - 56 (triangles) and 57 - 5g (crosses), as functions of different 
values of the couplings A', A" and A'" as indicated in the abscissa, 
with the coupling A fixed at the value A = 0.05. The interactions are 
of the XX type with J\ = J\ = K* = K } ( = J, Vj, k e B and V£ 6 5 , 
and the same results hold for ferromagnetic and antiferromagnetic 
couplings. Bottom (lower panel): the corresponding energy gap. 



collection of perfect Bell singlets. However, in the multi-pair 
case the energy gap decays rapidly as the number of pairs in- 
creases, thus making this configuration extremely sensitive to 
thermal fluctuations. 



2. Multipartite surface entanglement 

We have seen that a structure of many entangled pairs, ap- 
proaching a dimerized configuration of many Bell singlets oc- 
curs when the coupling parameters Aj between the different 
surface spins and the bulk of the network are significantly dif- 
ferent from each other. Let us now consider the opposite sit- 
uation, in which essentially all the surface spins interact with 
the bulk spins with the same weak coupling: Aj = A <SC 1, 
VjeS. 

In this case, as discussed in Sec. [TTJ the effective Hamilto- 
nian for the surface spins dynamics is that of a fully connected 
system, in the class of the Lipkin-Meshkov-Glick model Isill . 
The symmetries of the geometry of the effective model imply 
that the ground-state correlations have to be shared between 
all the spins, hence resulting in a long-range pattern of multi- 
partite entanglement. 

In order to gain insight into the structure of multipartite 
entanglement in the reduced ground-state density matrix of 
the surface spins (which is a multipartite mixed state), we 



FIG. 8: Top: quantum spin networks with square (left) and cubic 
(right) configurations, each with four surface spins and with ferro- 
magnetic interactions of the XX type (J x . k = J y . k = K* = K y ( = J, 
Vy, k 6 B and W 6 5). Bottom: fidelity between the multipartite 
entangled state in Eq. {8} and the reduced ground-state density ma- 
trix of the surface spins as a function of the surface-to-bulk coupling 
strength A. Inset: behavior of the corresponding energy gap as a 
function of A. 

have to bypass the problem of the lack of well defined and 
faithful measures of multipartite entanglement in multipartite 
mixed states. We will proceed first by introducing a reason- 
able ansatz (see Eq. ([8j for the pure ground state of the net- 
work's surface in the limit of small/vanishing coupling to the 
bulk, and we will discuss the structure of the associated resid- 
ual tangle, which is a well defined and faithful measure of 
multipartite entanglement for pure states 115 5fl . We will then 
show how the reduced ground state of the network's surface 
indeed approaches the proposed ansatz state as the coupling 
A of the surface to the bulk is progressively reduced. If for 
instance we consider a network with XX ferromagnetic inter- 
actions with a surface composed by four spins as illustrated in 
Fig. [8] then, following the discussion in Sec. |IlJ the effective 
surface dynamics is described by a fully connected effective 
XX Hamiltonian for the four surface spins; therefore the natu- 
ral pure multipartite ansatz for the network's surface is given 
by the ground state of the effective, fully connected surface 
Hamiltonian. For symmetric configurations, such as those il- 
lustrated in Fig. [8] each spin of the surface will effectively in- 
teract with all the other surface spins, with the same coupling 
strength. It is straightforward to verify that the ground state of 
the corresponding effective Hamiltonian, in the ferromagnetic 
case, is the following multipartite entangled state: 

|Z ) = -L (| TTU> + 1 UU> + 1 TUT> 
V6 

+1 mi> + 1 mt> + 1 utt» ■ (8) 
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By symmetry considerations, the ground state of the effec- 
tive 4-spin XX Hamiltonian must have net zero magnetization 
on the XY plane, and state |Zq), Eq. ||8), is the unique state 
of lowest energy with this property, in accordance with the 
fact that for gapless systems the ground state is non degen- 
erate. The normalized residual tangle, taking values in the 
interval [0, 1] lf55ll can be computed exactly for the state |Zo), 
Eq. ([8]), yielding 2/3 for each spin [56], a large finite value 
that demonstrates the existence of strong multipartite entan- 
glement between the surface spins. In the limit of very small 
values of the surface-to-bulk coupling A, state |Z ), Eq. ([8]), 
approximates closely the reduced ground sate po of the sur- 
face spins. This is shown in Fig. [8] where we report the fidelity 
F = Tr{po|Z())(Zo|} as a function of the parameter A. This be- 
havior demonstrates the presence of strong genuine multipar- 
tite entanglement among the surface spins in networks with 
highly symmetric interactions. 

IV. CONCLUSION AND OUTLOOK 

In the present work we have introduced and discussed the 
properties of surface entanglement, a generalization of the 
phenomenon of long-distance entanglement in quantum spin 
chains I2H - I24I1 to higher dimensional arrays of quantum spins. 
Surface entanglement is defined as the entanglement present 
in the reduced ground state of distant and non directly inter- 
acting spins belonging to the external boundary of a quantum 
spin network. The conditions for the occurrence of sizeable 
surface entanglement are that the surface spins be weakly cou- 
pled to the bulk spins and that the ground state of the bulk 
be non degenerate. We have observed that, typically, geo- 



metrically frustrated networks exhibit a weaker surface entan- 
glement compared to the corresponding geometrically unfms- 
trated ones. We have also discussed how the surface entan- 
glement is enhanced in networks with a modular structure, 
thereby exten ding the concept of one-dimensional modular 
entanglement 112511 to structures with modules of arbitrary ge- 
ometry and dimensionality. 

While long-distance entanglement allows to entangle only 
two distant spins at the ends of a quantum spin chain, sur- 
face entanglement may involve large numbers of spins result- 
ing in a great variety of entanglement patterns. Indeed, we 
showed that it is possible to entangle in parallel several spins 
pairs at the surface of the network or to create strong mul- 
tipartite entanglement between all the spins in the surface. 
These extended entanglement properties in the outer regions 
of a network with weak boundary-to-bulk links are promis- 
ing resources for routing quantum information among many 
distant nodes. Indeed, the ground-state surface entanglement 
that we have demonstrated is the result of effective interac- 
tions between the surface spins that are mediated by and have 
the same symmetries of the strongly correlated bulk. In this 
context, an interesting and relevant question is whether it may 
be possible to control the effective surface dynamics for quan- 
tum information and communication purposes, in particular 
for the efficient, parallel state transfer between many nodes in 
a network JUSSIH]. 
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